Generic susceptibilities of the half-filled Hubbard model in infinite dimensions 
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Around a metal-to- insulator transition driven by repulsive interaction (Mott transition) the single 
particle excitations and the collective excitations are equally important. Here we present results for 
the generic susceptibilities at zero temperature in the half-filled Hubbard model in infinite dimen- 
sions. Profiting from the high resolution of dynamic density-matrix renormalization at all energies, 
results for the charge, spin and Cooper-pair susceptibilities in the metallic and the insulating phase 
are computed. In the insulating phase, an almost saturated local magnetic moment appears. In the 
metallic phase a pronounced low-energy peak is found in the spin response. It is the precursor of 
the magnetic moment in the insulator. 

PACS numbers: 71. 27.-(-a,71.30.-fh,75.20.Hr, 71.28. -|-d 
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I. INTRODUCTION 

Strongly correlated systems persist to be a very inter- 
esting field of current research. In particular in the vicin- 
ity of a quantum phase transition the physics is very rich 
because the nature of the ground state and of the excita- 
tions changes. One prominent example of such a phase 
transition at zero temperature is the metal-to-insulator 
transition driven by an increasing local repulsive inter- 
action: the so-called Mott transition. For low values of 
the interaction the system is metallic because the elec- 
trons can still pass one another. For large values of the 
repulsive Coulomb interaction, there can be at most one 
electron per site. If the system is exactly half-filled, each 
site is occupied by one electron. No motion of electrons 
is possible because they are blocking one another. Hence 
the system is insulating. 

The simplest model describing the Mott transition is 
the Hubbard model. ^'^''^ The one-dimensional case can 
be solved analytically and has been studied intensively.^ 
It is governed by the particular phenomena of one- 
dimensional physics such as spin-charge separation. Yet 
this is not the generic physics occurring in higher di- 
mensions. Very much of our current understanding of 
the Mott transition in higher dimensions is based on the 
limit of infinite dimensions''^ which leads to the dynamic 
mean-field theory (DMFT) (Refs. 7 and 8) as important 
approximation scheme for real narrow-band compounds. 

The essential result of DMFT is that the transition 
between metal and paramagnetic insulator is marginally 
first order. ^'^'^^'^-'^'^^'^^ It is first order at finite temper- 
ature where a finite amount of spectral weight is redis- 
tributed at the transition Uc with Uca{T) < Uc{T) < 
Uc2{T). But at zero temperature only an infinitesimal 
amount of spectral weight is redistributed at Uc{T = 0) = 
Uc2{T = 0) (Ref. 13) so that the transition is continuous. 
The first order jump has just vanished. The insulator 
represents a metastable phase for U < Uc2 (Refs. 11 and 
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FIG. 1: (color online) Phase diagram of the Mott transition of 
a half-filled Hubbard model in infinite dimensions with semi- 
elliptic density-of-states at zero temperature as function of the 
interaction U . The critical value below which the insulator 
ceases to exist is Uc\ ~ 2.38±0.02D. The critical value above 
which the metal ceases to exist is Uc2 = 3.07 ±0. ID (Refs. 13 
and 14). But between Uci and Uc2 the metal has the lower 
energy so that the insulator is only metastable (Ref. 13). 



13) (see Fig. 1). 

Historically, there are two scenarios for the Mott tran- 
sition based on opposite limits. The Brinkmann-Rice 
scenario captures the essential point on the metallic 
side, namely the band- narrowing.^' The Hubbard sce- 
nario captures the two bands in the insulator (Hub- 
bard bands), which approach each other till they touch 
at the point where the insulator becomes unstable. 
The DMFT combines the strong points of both preced- 
ing scenarios. ^'^'""^"'"'^^'^'''^^'^^ Already the metallic solution 
displays Hubbard bands. The difference between the 
metallic and the insulating solution is found in the re- 
distribution of spectral weight at moderate energies only 
while the spectral densities at higher energies coincide. 
This coincidence is quantitative for U Uc2 as shown in 
Ref. 13. 

So the generic single-particle dynamics as encoded in 
the single-particle propagator is by now well-understood 
for the Mott transition. The dynamics of collective exci- 
tations and in particular the interplay between the single- 
particle modes and the collective modes is less well un- 
derstood and investigations are ongoing. ^^'^° 

A truly open issue is the physical origin of sharp fea- 
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turcs found at the inner band edges of the Hubbard bands 
just before the system switches from metallic to insu- 
lating behavior. These features were observed in quite 
a number of investigations but they were discussed as 
physical phenomenon only recently ^'^'^^ based on high- 
resolution dynamic density-matrix renormalization (D- 
DMRG). The features are confirmed by independent D- 
DMRG calculations^^ and by high-resolution numerical 
renormalization (NRG) (Ref. 22) but called into question 
by quantum Monte-Carlo extrapolations.^^ 

For the above reasons we have performed a thorough 
investigation of the susceptibilities in the half-filled Hub- 
bard model in infinite dimensions. The calculations start 
from the self-consistent solutions obtained by iterating 
the DMFT self-consistency cycle®'24.25 ^ith D-DMRG as 
impurity solver. ^'^'^'^ We benefit again from the good con- 
trol of the energy resolution achievable by D-DMRG for 
all energies. 

The susceptibilities provide valuable complementary 
information to the single-particle propagator. They ad- 
dress bosonic observables such as the local charge or spin 
or Cooper-pair density and their dynamics. So they give 
information about the corresponding collective modes. 
Moreover, the susceptibilities are experimentally rele- 
vant. The charge susceptibility corresponds to the po- 
larizability which determines the response seen in linear 
optics such as in infrared absorption. The spin suscep- 
tibility can be measured by inelastic neutron scattering. 
For these reasons, the susceptibilities ^'^^^'^s'^^ -^spm^ 
^pair g^j^g addressed in the present article. 

The paper is organized as follows. In Sec. II, we intro- 
duce the model. In Sec. Ill, the susceptibilities will be 
defined in detail. Analytic statements about them will 
be derived there. The numerical results in the insulat- 
ing and in the metallic phases will be given and their 
physical implications will be discussed. The conclusions 
summarize the article. 



II. MODEL AND METHOD 

We study the Hubbard model in DMFT at half-filling. 
The Hamiltonian reads 

^ = -i E 4<xS;. + ^E("^;T-l/2)(na-l/2). (1) 

Here c|.g. creates a fermion of spin a at site i while Ci-a 
annihilates such a fermion. The matrix element t labels 
the hopping of the fermions from site to site, i.e., their 
kinetic energy. The matrix clement U > labels the on- 
site Coulomb repulsion of two electrons on the same site. 
It is an effective parameter which takes the screening into 
account. 

We assume that the lattice on which it is defined is bi- 
partite so that the Hamiltonian (1) displays particle-hole 
symmetry. In our physical discussion we will assume that 
the system is translationally invariant so that the mo- 
mentum dependence of a propagator can be considered. 



The actual calculations, however, will be done for a non- 
interacting semi-elliptic density-of-states (DOS) 

poiiu) = {2/{^D'-)Wd^-lo^ (2) 

which is characteristic for the Bethe lattice with infinite 
branching ratio. The advantage of this approach is cal- 
culational simplicity and a finite support in contrast to 
the Gaussian tails on truly hypercubic lattices. ^'^ 

The limit of infinite dimensions d or equivalently of 
an infinite coordination number z exists if t is scaled 
like r/Vi. This hmit defines the DMFT.5'6.8,24,25 ^j^g 
self-energy E^j becomes local E^ — SijYin and equals 
the self-energy of a single-impurity Anderson model^^ 
which has the same skeleton diagrams in all orders. This 
means that the local dressed propagator Gu must be the 
same which defines the self-consistency condition of the 
DMFT. The ensuing simplification is that one only has 
to solve a zero-dimcnsional problem: an interacting site 
coupled to a bath. One way to represent this bath is 
as semi-infinite chain so that the problem is amenable 
to powerful one-dimensional tools, for instance dynamic 
DMRG^'^'2^ which ensures a good control of the energy 
resolution over all energies. For details, we refer the 
reader to Ref. 14. 



III. SUSCEPTIBILITIES 

In a translationally invariant system it is appropriate 
to consider the momentum-dependent x(q)- But the im- 
plications of the limit of infinite dimensions are easier 
seen in real space. In real space x depends only on the 
difference — Yj between the site i where the observable 
is measured and the site j where the field is applied. In 
the limit d — > cx) each fermionic propagator from i to j 
is scaled by a factor d"!!'"-'"/^ where || • || is the taxi cab 
or New York metric which counts the minimum number 
of hops required to get from i to j. In a diagrammatic 
description of the propagation of any bosonic, collective 
observable from i to j at least two fermionic propagators 
link i and j. This implies that such a susceptibility is 
suppressed by d""'"-'!!. Hence only local susceptibilities 
from i to i matter in infinite dimensions, at least in the 
absence of phase transitions. 

This conclusion is not quite the whole story because 
certain non-local contributions can add up. For instance 
there are 2c? next-neighbor contributions (||i — j|| = 1) 
to i so that they make a non-negligible contribution if 
they add up since 2d/d ~ 2. But they will not add up 
for a generic momentum q. The momentum q enters the 
susceptibilities only via 

1 

V{<i) = -^^^coaiqp) (3) 

p=i 

where qp is the component in direction p.^'^" In the 
limit d — > cx) almost all vectors q imply 77(q) = 
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since is gaussian distributed with finite vari- 

ance. Only particular values of measure zero, for instance 
q = (tt, TT, TT, . . . , 7r)t, imply a non-vanishing ry. Hence, 
the generic susceptibility is the one for rj = which cor- 
responds in real space to the local one. 

Of course, some important effects of finite-dimensional 
physics are not captured by the generic susceptibilities. 
For example an antifcrromagnctic instability or the in- 
stability to an incommensurate phase is indicated by the 
divergence of a susceptibility for some ?7 7^ 0. This must 
be kept in mind. On the other hand, however, the prop- 
agation of collective modes as they interact with single 
particles is given in c? = 00 by the generic susceptibilities. 
In any diagram for the proper self-energy which contains 
the propagation of a collective mode there is also a sum 
over its momentum. Hence the peculiar contributions 
with rj ^ do not matter here. It is the generic behav- 
ior of collective modes at 7/ = which is relevant for the 
interaction with single particles in d = 00. 

Note that this argument remains true even if the 
DMFT is not seen as the limit of infinite dimensions but 
more broadly as a consistent local approximation scheme. 
This is by now a very common view adopted in the de- 
scription of real compounds. Then one can discuss the 
whole momentum dependence of the susceptibility. But 
the behavior of the collective modes which enters implic- 
itly in the description of the single-particle dynamics re- 
mains the one given by the local susceptibilities. For this 
reason, we compute and discuss the local susceptibilities 
in the following. 

The local susceptibilities are easily accessible since 
they are identical to the local susceptibilities at the in- 
teraction site of the auxiliary single-impurity Anderson 
model. This is obvious if one thinks of the susceptibilites 
as being given by an expansion in terms of skeleton dia- 
grams to infinite order. 

In the sequel, we will present the imaginary parts of 
the susceptibilities since these parts provide direct infor- 
mation on the energies and spectral weights of collective 
excitations. 

One drawback of the susceptibilities is that their imag- 
inary parts are antisymmetric (odd) by construction 
Imx([j) = — Imx(— w). Hence the very interesting be- 
havior at zero and at very low energies is suppressed. 
For instance, let Q be an hermitean bosonic observable 
like the spin density. With |0) being the ground state 
and Eq its energy, the matrix element of the resolvent 

R{u;) := {0\QiLj- (n- Eo))-^Q\0) (4) 

can have a 6 peak at = in its imaginay part. This 
would constitute an important piece of information on 
the system. But the susceptibility x constructed from R 
according to 

x{uj)^-R{oj + iO+)-R{-uJ-iO+) (5) 

would not display this 5 peak because it cancels on the 
right hand side of (5). In order not to lose the informa- 
tion at zero and at very low frequencies we will display 



results for 

lmX+{Lu) = lmRiuj) (6) 

which is the contribution for non- negative frequencies. 
The contribution for negative frequencies is implied by 
antisymmetry, i.e., Imx-(w) = ~lmx+{—to)- 

An important tool in understanding spectral densities 
are sum rules. From the relation (6), the definition (4), 
and the Hilbcrt representation of R{uj) it is obvious that 
the total spectral weight takes the value 

/•OO 

/ x+i^)di^ = TT lim uoR[uj) (7a) 
= ^(0|Q2|o). (7b) 

So knowing the ground state expectation value of 
helps to understand general trends in spectral weights 
and it provides an important check for numerical calcu- 
lations 

A last important point concerns the computation of 
spectral densities like Imx+(ct') by D-DMRG. The dy- 
namic DMRG calculates a correction vector |cv) besides 
the ground state |0) and the state obtained from the ap- 
plication of Q, QjO).'^^'"^^ This correction vector reads 

\cy)^[lo + i5-{H-Eo))-^Q\Q). (8) 

Obviously, (0|Q|cv) yields -R(tJ-f-i(5). The computation of 
the correction vector requires a numerically demanding 
matrix inversion. If there is an inaccuracy e in the cor- 
rection vector |cv) the imaginary part of R{uj + i5) can be 
obtained from a variational functional with a decreased 
inaccuracy of the order of |ep, see Ref. 33. 

In any case, the numerical approaches require the 
imaginary frequency 5 in (8) to be finite. This implies a 
certain broadening of the actual spectral density. In or- 
der to retrieve the unbroadened spectral density Imx+ (w) 
we employ the non-linear least-bias (LB) deconvolution 
technique.'^'* This approach yields always non-negative 
results as is to be expected for Imx-|.(a;). For details of 
the approach wc refer the reader to Rcf. 34. Any devia- 
tions from the procedure described therein will be given 
below where applicable. 

A. Cooper Pair Susceptibility 

Here we consider the local observable 

which creates or annihilates a Cooper pair on site i. Of 
course, one would expect that the corresponding suscep- 
tibility x^'^"{'^) is strongly suppressed in a Hubbard band 
with repulsive interaction. Yet it can contain interesting 
features, for instance at higher energy. 

But it is not necessary to compute x^^" (1^) separately. 
Indeed, there is an underlying symmetry which links 
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QP''" to Q'='''^'■8^ see Eq. (12) below. As a result the 
pair susceptibility is identical to the charge susceptibil- 
ity. So we will not discuss x^^"i^) here but refer to the 
next subsection where x'^^'^'^^'^i^) is investigated. 

The symmetry becomes apparent under the transfor- 
mation cj. 7^ according to 

7I := c| cos (/3 — C| sin (/? (10a) 

7I :~ c| cos(y5 -f- C| sin(p. (10b) 

On a bipartite lattice this transformation is performed 
on all even sites with the angle f and on all odd sites 
with the angle —tp. Then the hopping terms in (1) are 
left invariant, independent of the value of ip. The same 
is true for the onsite interaction as given by the term 
proportional to C/ in (1). So the total Hamiltonian (1) 
remains invariant under (10). 

The interesting relation is the one for the observables 

Qchargc ^ gchargc cos(2(/j) + Q^^" sin(2(^) (11) 

where we use 

gcharge anticipation of Eq. (12). The 
subscript c refers to the expression in terms of the origi- 
nal fermions c and c' while the subscript ^ refers to the 
expression in terms of the transformed fermions 7 and 
7^. Here the value of the angle of rotation tp matters. 
For = 7r/4 we switch from the charge observable to the 
pair observable. Hence the corresponding local suscepti- 
bilities are indeed the same. No additional numerics is 
needed in the bipartite half-filled case. 

The above symmetry transformation has not gone un- 
noticed. It is one of the transformations at the basis of 
the SO (5) theory which is comprehensively reviewed in 
Ref. 35. One further important conclusion is that charge 
and superconducting order are degenerate as far as they 
stem from a bipartite Hubbard model at half-filling, for 
instance with negative U. 

B. Charge Susceptibility 

Here we consider the local observable 

Q^^^^^^^,^^^+^^^,^^^_^ (12) 

which measures the charge fluctuations around half- 
filling, i.e., the deviation of the total fermion number 
per site from 1. 

The general sum rule (7) requires the expectation value 
(0|(5^|0) which amounts for Q^hargc twice the dou- 

ble occupancy value. First we note 

(Qchargc)2^2nTni-(nT+ni) + l (13) 
where ha — cta^i-a- This imphes at half-filhng 

POO 

/ xt"'^"HdoJ = 271 {0\h^hi\0) . (14) 
"'0 

So the static quantity to be known for the sum rule is 
the double occupancy (see Refs. 36 and 13). 
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FIG. 2: (color online) Positive imaginary part x'f^'^^'^ of the 
local charge susceptibility in the insulating phase as function 
of frequency for various values of the interaction U. The arrow 
points in the direction of increasing interaction. 



1. Insulator 

In Fig. 2 a series of numerical results is shown for the 
positive imaginary part ^'^'^''so ^-j-^^ local charge sus- 
ceptibility in the insulating phase. In the D-DMRG we 
kept m = 128 states in the truncated DMRG basis. The 
mesh of frequencies is given by the interval Auj = 0.05Z? 
and the imaginary broadening was S = O.ll?. The LB 
dcconvolution'^^ was perforemd with a tolerance constant 
of 1/Alb = 1/100. The curves are not perfectly smooth 
but display some wiggles. This is due to the deconvolu- 
tion procedure employed. 

Physically, no special features are discernible. But 
two trends are clearly visible. First, the susceptibility is 
more and more suppressed as the interaction is increased. 
This results in an overall reduction of the area under the 
curves. It can be quantified by the sum rule (14). So it 
is natural that the spectral weight of the charge response 
decreases on increasing U because the latter suppresses 
the double occupancy (0|n|nj0) more and more.^'^^'^^ We 
have checked this sum rule numerically and found it to 
be fulfilled to within a relative error of 1.8% on the de- 
convolved data in the interval uj G [0, 5.9D]. 

Second, the spectral weight is shifted to higher and 
higher frequencies on increasing interaction. This is seen 
in two features. One is the peak position which is shifted. 
Its shift corresponds to the shift of spectral weight in the 
single-particle propagators. ^^'^^'•^^ These shifts reflect the 
simple fact that the energy difference between the lower 
and the upper Hubbard band is given by about U. Hence 
it increases linearly with U. 

The other feature is the onset of finite spectral den- 
sity which increases also with U. Due to the LB 
dcconvolution'^^ there is no region where the spectral den- 
sity is strictly zero. But we have checked that the sus- 
ceptibility data is perfectly consistent with the natural 
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FIG. 3: (color online) Positive imaginary part t^.^^'^sc ^j^^ 
local charge susceptibility as function of frequency for various 
values of the interaction U in the metallic phase. The arrow 
points in the direction of increasing interaction. Note the 
different scale of the response compared to the response in 
the insulator in Fig. 2. 

assumption that the onset of the Imx'^^'^^'^ takes place 
at 2A vifhere A is the single-particle gap (for data, see 
13 and 14). For this check we analyzed the susceptibility 
data by fitting a quadratic onset oc (cj — Wonsct)^ plus 
higher order corrections to its continuum. The quadratic 
onset is to be expected from the square root onset of 
the single-particle bands^'^''^^ which is convolved with it- 
self in the standard particle-hole bubble. Because the 
single-particle gap rises upon increasing U the collective 
response is shifted to higher and higher energies. 

The onset at 2A reflects the fact that the collective 
mode is made from a particle and a hole. In the particle- 
hole symmetric case considered here both have the same 
gap A so that any collective continuum starts only at 
2A. A lower onset, i.e., at lower energy, could arise only 
if binding transferred spectral weight to the frequency 
interval below 2A. No indication for such a binding is 
found here. 



2. Metal 

In Fig. 3 a series of numerical results is shown for 
the positive imaginary part x+^'^^° of the local charge 
susceptibility in the metallic phase. In the D-DMRG 
we kept between m = 128 to m = 256 states in the 
truncated DMRG basis; the frequency mesh is given by 
Aw = 0.05D and the imaginary broadening by S — O.ID. 
The LB deconvolution"^^ was performed with tolerance 
constants 1/Alb = 1/10 and 1/Alb = 1/100. The curves 
are not perfectly smooth but display some minor wiggles. 
This is due to the deconvolution procedure employed. 

As in the insulating regime we find the trend that in- 
creasing interaction suppresses the charge response. This 



is expected because it is related to the same sum rule (14) 
which holds independent of the phase under study. Fur- 
thermore, the spectral response is shifted to higher and 
higher energy. Again this general trend can be related to 
the same trend in the single-particle propagators. ^'^'^^ 

Interestingly, the charge response in the metallic phase 
displays much more structure than in the insulating 
phase. For low values of U we find a linear increase 
with frequency to for not too high values of ^ 0.61?. 
This is the expected behavior for a Fermi liquid. Its 
slope becomes smaller and smaller as the quasiparticles 
become heavier and heavier. From U ~ 1.5D onwards, 
most of the charge response lies in an intermediate range 
D ^ u! ^ 2D. There is still some spectral weight at 
lower frequencies but it is decreasing rapidly. An ad- 
ditional shoulder situated between 0.61? and ID occurs 
above U = 1.81?. We will discuss this feature in detail 
below. Above U ~ 2.2D there is a third rather flat hump 
discernible centered around uj — 3D. 

Qualitatively, the three regions of charge response can 
be understood from the single-particle response. The 
single-particle response, see for instance Figs. 12 and 13 
in Ref. 14, is mainly characterized by the heavy quasi- 
particle in the narrow central peak and by the broad 
emerging Hubbard bands of signiflcant weight which are 
centered around u w 1.51?. Assuming that the col- 
lective response is roughly given by a single diagram- 
matic particle-hole bubble (or by an analytic function of 
this bubble as in the random phase approximation or in 
more sophisticated approaches such as the local moment 
approach), '^^''^^ we simply have to convolve the single- 
particle response at positive frequencies with the one at 
negative frequencies. 

The response at low frequencies u ^ 0.81? stems from 
the convolution of the central peak of heavy quasiparti- 
cles with itself. It dominates at low values of U, but on 
increasing U it decreases in weight like Z'^ where Z is 
the quasiparticle weight vanishing linearly for U — * Uc2- 
The quasiparticle weight Z measures the spectral weight 
in the central peak. 

The response at intermediate frequencies 0.81? ~ w ^ 
2.21? stems from the convolution of the central peak of 
heavy quasiparticles with one of the Hubbard bands. 
Hence it is higher in frequency, because the Hubbard 
band is, and its weight decreases only linearly in Z 

The response at higher frequencies results from the 
convolution of the lower and the upper Hubbard band. 
Hence it is located around twice their energy, i.e., around 
3Z?. This contribution is not suppressed by Z so that it 
survives in the limit U —^ Uc2- This is consistent with 
our data for U = 2.8D shown in Fig. 3. Unfortunately, 
no reliable data even closer to the critical value Uc2 could 
be obtained. 

Note that only the last contribution at the higher fre- 
quency has an analog in the response in the insulator 
since there only the Hubbard bands exist. Indeed, the 
insulating response describes the high frequency metallic 
response very well as is illustrated in Fig. 4. The insulat- 
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FIG. 4: (color online) Positive imaginary part x+^^'^^'^ of the 
local charge susceptibility as function of frequency for inter- 
action U = 2.8D in the metallic and insulating phase. 



^ 0.4 



' 1 ' 1 ' 1 


' 1 


1 1 1 1 1 1 1 

metal ~. 




\ 
\ 
\ 


U = 2.0 D 




U = 2.4 D - 




M 


U = 2.6 D '_ 






U = 2.8 D 












K ^ ~ 






\ ^ ^ 




1 





0.5 1 1.5 2 2.5 3 3.5 4 

co/D 

FIG. 5: (color online) Positive imaginary part t^-^'^'sc ^j^^ 
local charge susceptibility as function of frequency for values 
of the interaction U in the metallic phase, but close to the 
transition to the insulator. The main contribution is extrapo- 
lated as if there were no additional shoulder (thin solid lines). 
Then the shaded area is attributed to spectral weight of the 
shoulder. 



ing and the metallic curve agree very well above w « 3D. 
The sizable differences below this frequency are remark- 
able in view of the shift of fairly little spectral weight 
between the metallic and the insulating single-particle 
solution. For experiment, for instance infrared absorp- 
tion. Fig. 4 provides valuable information how different a 
metallic and an insulating system can look even though 
only tiny parameter changes are made. In the present 
case, even no parameters are changed, but only different 
hysteresis branches are considered. 

Let us come back to the shoulder seen between uj w 
0.6D and lo w \D for U ^ 1.8D. Its position corre- 
sponds very precisely to the frequency where the sharp 




QP weight Z 

FIG. 6: (color online) Spectral weight W attributed to the 
shoulder in the metallic charge response as function of the 
quasiparticle weight Z. The solid line is a quadratic fit with 
c = 0.758. 



feature at the inner band edges has been found, sec Fig. 
2 in Ref. 13 and Figs. 12 and 13 in Ref. 14. Hence it 
is to be expected that there is a relation between both 
features. In view of the hypothesis that the sharp feature 
is caused by a resonance made from a heavy quasiparti- 
cle and a collective mode (see Refs. 13 and 14), it would 
be appealing to interprete the shoulder in Fig. 3 as the 
cause for the sharp feature in the single-particle spectral 
density. To pursue this idea further we plot in Fig. 5 the 
metallic charge response close to the transition to the 
insulator. 

In the curves shown in Fig. 5 the shoulder is clearly vis- 
ible. We extrapolate the main peak on which the shoul- 
der sits smoothly. This is done by determining a fre- 
quency interval [6i , 62] below the shoulder and a second 
one [63,64] above it where the shoulder is not present. 
These intervals are found from analyzing minima and 
points of inflections of the original curve, for example for 
U = 2 AD we took [0.56£i, 0.646£»] and [1.038Z3, 2.1751?]. 
Then the data within these two intervals is interpolated 
by a 10*'' order polynomial. This is taken to describe 
the continuum without the shoulder. The weight of the 
shoulder (shaded area in Fig. 5) is given by integrating 
the difference between the original curve with shoulder 
and the lO"' order polynomial in the interval [62, 63]. The 
resulting weights are well-defined within 4 • 10~^. 

This procedure is applicable for the results obtained 
for UjD e {2.0, 2.2, 2.3 . . .2.7, 2.8}. The resuhing values 
are depicted in Fig. 6 as function of the quasiparticle 
weight Z . The quadratic fit agrees very well with the 
data except for the last point resulting from U ~ 2.QD. 
For such a fairly low value of U the separation of the 
shoulder from its background is not possible reliably. 

Clearly, the shoulder weight W depends quadratically 
on Z. We recall that the spectral weight S of the sharp 
feature at the inner band edges in the single-particle 
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propagator depends linearly on Z: S (x Z as found 
previously. These facts are incompatible with the 
sharp feature S being the result of the shoulder W. It 
would require that right at the transition to the insula- 
tor the shoulder induces the sharp feature although the 
weight of the shoulder is infinitely smaller than the weight 
in the sharp peak. 

But the other way around the quadratic behavior in 
Fig. 6 finds its natural explanation. The shoulder results 
from the convolution of the central quasiparticle peak of 
weight Z with the sharp feature with S oc Z. Hence 
W (X Z^ ensues as found. 

Also the position in frequency is explained in this way. 
Since the central peak is located at zero frequency the 
shoulder as result of the convolution with the sharp fea- 
ture is located at the frequency where the sharp feature 
is found. 

Summarizing these findings, we conclude that the 
shoulder in the charge response can be understood as a 
consequence of the sharp feature at the inner band edges 
of the metallic single-particle spectral density. It is not 
its cause. While it is satisfying to have explained the 
origin of the shoulders in Figs. 5 and 6 we state that the 
physical origin of the sharp feature described in Refs. 13 
and 14 is still unresolved. 



C. Spin Susceptibility 

Here we consider the local observable 

Q^'^'" = 4T^.;T-4iC.a (15) 

which measures the spin fluctuations around zero mag- 
netization in z direction. 

The general sum rule (7) requires the expectation value 
(0|Q^|0) which amounts for Q^^'" up to an expression 
which contains again the double occupancy. First we 
note 

(Q^P'")2 ^n^+h^- 2ntni (16) 

where we used that the fermionic occupation number h 
is equal to its square. This implies at half-filling 

/ xf"Hrf^ = 7r(l-2(0|ritrix|0)). (17) 
Jo 

Note that this expression stays finite in the limit of van- 
ishing double occupancy as it occurs for U ^ oo. 

1. Insulator 

In Fig. 7 a series of numerical results is shown for the 
positive imaginary part of the local spin suscepti- 

bility in the insulating phase. In the D-DMRG we kept 
m — 128 states in the truncated DMRG basis the mesh 
is given by the frequency interval Alu — 0.05D, and the 
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FIG. 7: (color online) Positive imaginary part x+"" of the 
local spin susceptibility in the insulating phase as function 
of frequency for various values of the interaction f/. The ar- 
row points in the direction of increasing interaction. The 
Lorentzian at zero frequency represents the example at f/ = 
4_D for the broadened 5 peak occurring in the insulator. 

imaginary broadening is 5 = O.ID. In the LB deconvolu- 
tion, the tolerance constant Alb = 1/10 is used.'^"' The 
curves are not perfectly smooth but display some very 
small wiggles due to the LB deconvolution. The attempt 
to deconvolve the numerical DMRG data as a completely 
continuous spectral density leads to a very large and very 
narrow peak at low frequency (not shown) . The continua 
beside this dominating term cannot be resolved reliably. 
But it turns out that the ansatz 

ilmxf" = A<5(c.)+p™-*(^) (18) 

TT 

works extremely well for deconvolution, see Fig. 7. Here 
p'=°"'(a;) stands for the continuous spectral density which 
is retrieved via the LB deconvolution. The weight A of 
the b peak results from the non-linear set of equations 
defining the Lagrange multipliers appearing in the LB 
ansatz.'^'' 

Why does a zero frequency 8 function make sense phys- 
ically in the local spin response of a paramagnetic insu- 
lator in infinite dimensions? The exchange coupling J 
in a Heisenberg model derived from a Hubbard model in 
the insulating regime reads J = At^ jU Hence scaling 
t = t*/y/z implies J <x 1/z and the exchange coupling 
does not contribute unless all the nearest neighbors con- 
tribute on average the same non- vanishing amount. This 
implies that a static mean-field treatment of the Heisen- 
berg antiferromagnet around the Ising limit becomes ex- 
act in infinite dimensions. There are no short-range 
spin-spin correlations. Hence each spin feels only the 
local field /imf = —zJm generated by the average mag- 
netization m of its z neighbors. But in the paramagnetic 
phase which we consider here the average magnetization 
is zero: m = 0. Hence there is no field and concomitantly 
there is no preferred direction of the local spin. This im- 
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FIG. 8: (color online) The total weight (circles) and four times 
the square of the local Sz component (squares) of the spin 
response x+""(t^) in the insulating (filled symbols) and in the 
metallic phase (open symbols). In the insulator the local Sz 
component is taken from the weight A of the 5 peak as in Eq. 
(18). In the metal it is taken to be the weight of the dominant 
low energy peak, see Figs. 9 and 10. The dark lines interpolate 
the data in the insulator while the light lines interpolate and 
extrapolate the data in the metal. 

plies that the spin response is the one of a free spin, this 
means, a b function at zero frequency. It signals that a 
local spin flip docs not cost any energy. 

Besides understanding this physics on the level of the 
infinite dimensional Hubbard model it can be under- 
stood on the level of the effective impurity model. In- 
deed, a gapped impurity model at particle-hole symme- 
try is found to be always in the so-called local moment 
regime, ^^'^^'^^ where a free local moment is formed. 

The weight A of the i5 function represents four times 
the square of the local S^, component in our normaliza- 
tion. It may not be confused with a static magnetiza- 
tion which is absent in the paramagnetic phase consid- 
ered. For \J oo, A takes the value unity since the 
spin 5' = 1/2 is fully localized. At any finite interaction 
t/ < oo, A is reduced since charge fluctuations renormal- 
ize its value downward. The spins are not fully localized 
but smeared out to some extent to adjacent sites due to 
their virtual excursions. This behavior is analyzed quan- 
titatively in Fig. 8 where the filled squares denote the 
values of A. Further discussion is presented below. 

Besides the b peak a continuous contribution of low 
weight persists. It correponds to the charge fluctuations 
which take some weight away from the dominant local 
spin response at zero frequency. Indeed, the continua arc 
very similar, though not identical, to the ones found in 



the charge response in the insulating regime, see Fig. 2. 
To the accuracy that our numerical deconvolved data al- 
lows the onset of the continuous spectral spin response 
is again at 2 A, i.e., twice the single-particle gap in the 
insulating regime. Hence the origin of the continuous 
spin response is essentially the excitation of particle-hole 
pairs. This is qualitatively similar to what one expects 
from the diagrammatic result in random-phase approxi- 
mation where the response is a function of the particle- 
hole bubble. 

The fact that most of the spectral weight in the spin 
response is found at zero frequency and not in the con- 
tinuum starting at 2A shows impressively that a binding 
phenomenon occurs. The spin response at zero frequency 
can be viewed as the signature of a bound state of a 
particle-hole pair.^* Only by binding one can understand 
how spectral weight can be transferred to energies lower 
than the sum of the energies of the constituent states. We 
think that such shift of spectral weight due to binding is 
not taken into account by the sophisticated argument on 
the spectral density close to the Mott transition.*"^ This 
argument excluded the continuous scenario at zero tem- 
perature which is supported by most other analytical and 
numerical evidence''^^ (see Sec. I for a sketch of this sce- 
nario) . 

The total weight of the spin response as displayed in 
Fig. 8 just reflects the behavior of the double occupancy 
according to the sum rule (17). Hence it approaches unity 
for [/ — > oo but it does not become very small ovlV ^ Ud 
either. 

Much more interesting is the behavior of the square of 
the local Sz component as quantified by A in Eq. (18). 
The square root of this expression can be identified with 
the local magnetic moment. Clearly, A ~ 1 a.t U ^ oo 
is the starting point. But it is remarkable that A w 0.94 
has hardly decreased for U = Uc2 where the insulator is 
no longer the ground state. In physical terms this means 
that a Mott insulator is governed by very well localized 
spins as long as it exists. Hardly any renormalization 
due to charge fluctuations takes place. From a theoreti- 
cal point of view this can be explained by the significant 
charge gap A w 0.45D at 11^2 (Ref. 13) which acts as 
an infrared cutoff limiting the influence of charge fluctu- 
ations. This can be easily understood by the renormal- 
ization flow of the impurity model where the insulator 
corresponds to the local moment fixed point. *^'*^ 

Even more remarkable is that the local magnetic mo- 
ment is not much lower at Uci either. Below this in- 
teraction the insulator is not longer locally stable. Even 
there A is still larger than about 0.87 although the charge 
gap has become zero and the lower and the upper Hub- 
bard band arc touching each other, sec for instance Fig. 
2 in Ref. 13. But it is obvious from the touching Hub- 
bard bands that the DOS p{ui) at the Fermi level, i.e., at 
w = 0, is still zero. We conclude that no hard infrared 
cutoff is needed and that the fact that the insulator at 
Uci displays a semi-metallic DOS with liuY^^^^j p{^) = 
is sufficient to bound the influence of the charge fluctu- 
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ations. So the magnetic moment is not rcnormalizcd to 
zero and the fixed point of the renormahzation of the 
corresponding impurity model is still the local moment 
fixed point. 

The main goal of this paper is the comprehensive 
analysis of the susceptibility around a Mott transition 
in infinite dimensions. But in view of the remarkable 
findings at d = oo it is in order to speculate how these 
findings change on passing to finite dimensional systems. 

The main difference is that any finite dimensional sys- 
tem would show at least short-range magnetic correla- 
tions. Hence the magnetic response would not be gov- 
erned by a i5 peak at zero frequency as in Fig. 7. If the 
resulting antiferromagnetic system is sufficiently strongly 
frustrated and/or sufficiently low-dimensional so that the 
magnetic fluctuations are strong enough to prevent mag- 
netic long-range order the system would be paramag- 
netic displaying a magnetic gap. Generically, the mag- 
netic excitations would be triplons^^ with some disper- 
sion. Hence the local spin response would show the sum 
of triplon contributions from all the wave vectors in the 
Brillouin zone. A broad feature at finite frequencies in a 
frequency range given by the magnetic exchange J would 
be seen in Imx^™(a;). A sharp mode at finite, but low 

frequency would be discernible in Imx+""(ci;, q) at a given 
wave vector q. The sum of the weights in these sharp 
modes over the Brillouin zone constitutes the finite di- 
mensional analog of the weight A in our infinite dimen- 
sional analysis. We expect other features to be qualita- 
tively very similar to the above findings. For instance the 
local magnetic moment in any insulating state should be 
very little renormalized due to charge fluctuations. 



2. Metal 



In Fig. 9 a series of numerical results is shown for the 
positive imaginary part x+"" of tl^c local spin suscepti- 
bility in the metallic phase for not too large values of the 
interaction. The curves for larger values of U are plotted 
in Fig. 10 In the D-DMRG we kept between m = 128 and 
m = 256 states in the truncated DMRG basis, the fre- 
quency mesh is given by the intervals Aw = 0.025Z? and 
Aw = 0.05-D, and the imaginary broadening is chosen be- 
tween 5 = 0.05D and 5 — O.ID. The LB dcconvolution'^^ 
is done with the tolerance constant 1/Alb = 1/10. 

As a first check of our data we compute the static spin 
spin susceptibility x''p™(0) via the Kramers-Kronig rela- 
tion 
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X^P'"(0) = - 

TT 



Imx+"(w) 



Imx+'"(-w) 



dui. (19) 



The negative term in the numerator occurs only be- 
cause the LB deconvolution tends to produce spurious 
minor contributions at negative frequencies; otherwise 
Imx+™(w < 0) is strictly zero at zero temperature. The 
results are compared to previous results in Fig. 11, see 
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FIG. 9: (color online) Positive imaginary part of the 

local spin susceptibility in the metallic phase as function of 
frequency for various small values of the interaction I] . The 
arrow points in the direction of increasing interaction. Note 
the quickly increasing peak at low frequencies and the quickly 
vanishing spectral weight at higher frequencies. 




FIG. 10: (color online) Positive imaginary part x+"" of the 
local spin susceptibility in the metallic phase as function of 
frequency for various values of the interaction U close to the 
transition to the insulator. The arrow points in the direction 
of increasing interaction. The peak at low frequencies prevails 
completely, note the scale on the y axis. 



Fig. 44 in Rcf. 8, which were obtained by exact diagonal- 
ization. The agreement is very good at low U deteriorat- 
ing for larger values of U . We attribute the discrepancy 
at larger values of U to effects of finite size and of finite 
temperature in the exact diagonalization approach. But 
for small and intermediate values of the interaction the 
comparison underlines the validity of our results. 

The curves are dominated by a prominent peak at low 
frequencies. In some curves, in particular between IJ = 
\D and 2_D, there appears to be a shoulder to this peak. 
Since this feature occurs already for moderate values of \J 
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FIG. 11: (color online) Static susceptibility x°^'"(0) from our 
dynamic data via (19) (black crosses) compared to data ob- 
tained by exact diagonalization at very small temperature 
(red circles), adapted from Fig. 44 in Ref. 8. The black 
cross at [/ = corresponds to the analytic result x'''""(0) ~ 
16/(37rD). Lines are guides to the eye only. 
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FIG. 12: (color online) Positive imaginary part x+"" of the lo- 
cal spin susceptibility as function of frequency for interaction 
U = 2.SD in the metallic and insulating phase. 



where the curves arc still fairly smooth wc arc confident 
that this shoulder is in fact a real physical feature. But 
due to the dcconvolution procedure involved we cannot at 
present rule out completely that it is a spurious numerical 
effect. In the following we refrain from discussing the 
shape of the low-frequency peak. 

At higher frequencies around lu « 2D in Fig. 9 a very 
broad peak of low amplitude can be seen. Qualitatively, 
this continuum results from the convolution of the lower 
and the upper Hubbard band or more complicated de- 
scriptions in terms of excitations from the lower and the 
upper Hubbard band.'^^'^® It corresponds to the contin- 
uum found in the insulating phase in Fig. 7. 

The insulating and the metallic response at higher 



frequencies are alike, see Fig. 12, because their single- 
particle spectral densities are identical for higher frequen- 
cies for U Uc2, see Fig. 2 in Ref. 13. This is analo- 
gous to what we have discussed in the charge response 
in Fig. 4. In both cases, the insulating and the metallic 
responses coincide for uj ^ 3D. The slightly more wiggly 
metallic spin response in Fig. 12 results from the difh- 
culty to deconvolve the relatively small continuum close 
to the dominating low-frequency peak. 

Since the low-frequency peak is clearly the dominat- 
ing feature its physical significance has to be elucidated. 
Clearly, it is shifted towards zero frequency on U ^ Uc2 
while becoming narrower and narrower. So it is to be 
expected that it represents the precursor in the metal of 
the 5 peak in the insulating regime. In order to support 
this claim we analyze the sum rule (17) and the weight in 
the low-frequency peak. Both sets of data are depicted 
for the metal in Fig. 8 by the open symbols. The sum 
rule is again well fulfilled to within a absolute relative 
error of 0.5%. The weight in the low-frequency peak is 
integrated till the first minimum in the spectral density 
on the right hand side of the peak is reached. First, we 
note that the total spectral weight of the metal equals 
the one of the insulator for U Uc2 as far as the ex- 
trapolation can be trusted. This is expected since the 
transition occurs precisely where the double occupancies 
become equal. ^^'^^ 

Second, wc note that also the weight of the dominant 
low-frequency peak in the metal approaches the weight 
of the 5 peak in the insulator to very good accuracy. 
This clearly corroborates our hypothesis that the metal- 
lic dominant low-frequency peak is the precursor of the 
bound state at zero frequency in the insulator. Naturally, 
there is no sharp bound state in the metal because such 
a state can decay into particle-hole states made from the 
heavy quasiparticles and -holes which are still present in 
the metal. Hence no bound state but a resonance oc- 
curs. The width of this resonance can be understood as 
Landau damping. 

Recall that the low-energy spin resonance in the gap- 
less case becomes the zero-frequency mode of the cor- 
responding gapped case in single impurity Anderson 
models. "^^ So our analysis is well-founded also on the level 
of the impurity models. 

In the insulating regime we have tentatively carried 
our infinite dimensional results over to finite dimensions. 
We have speculated that the 5 peak in Fig. 7 becomes a 
dispersive magnon, if magnetic long-range order exists, 
or a dispersive triplon if not. Both, magnon or triplon, 
are bound states of particle-hole pairs from the electronic 
point of view. So we expect that the emergent magnetic 
resonance found here becomes in finite dimensions a dis- 
persive resonance which is the precursor of a perfectly 
sharp magnetic excitation. In literature, the term 'para- 
magnon' is used for such precursive resonances. If no 
magnetic order is to be expected the term 'paratriplon' 
would be more appropriate. 

We recall that the existence of procursive magnetic ex- 
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citations and their interaction with the single-particle 
excitations is very important for the understanding of 
kinks in the electronic dispersions^^'^" and possibly also 
for Cooper pairing in strongly interacting systems. 

IV. CONCLUSIONS 

In this paper, we investigated the zero temperature 
Mott transition as function of the interaction in a generic 
model, namely the half-filled Hubbard model in infinite 
dimensions. We focused on the susceptibilities which are 
of theoretical and experimental relevance. Thereby, com- 
plementary information to the existing investigations of 
the single-particle dynamics is provided. 

We showed that in infinite dimensions the generic sus- 
ceptibilities are the local ones. Locally, only three types 
of bosonic observables exist: the charge, the spin and 
the Cooper pairing operator. So we discussed the cor- 
responding susceptibilites T^>^hargc^ X^^'", and x^^" ■ By 
an intricate symmetry, x^^" is found to be identical to 

T^chargc 
A 

For the charge susceptibility x'^^'^'^^'^ in the insulating 
phase we found a strong suppresion on increasing repul- 
sive interaction. The spectral density sets in at 2A, i.e., 
at twice the single-particle gap. No binding phenomenon 
occurs; the spectral line is rather featureless. 

In the metallic phase, three ranges in frequency can be 
distinguished. The first results from a heavy quasipar- 
ticle and a heavy quasiholc, the second from one heavy 
excitation and one in one of the two Hubbard bands, 
and the third consists of a particle in the upper and a 
hole in the lower Hubbard band. On approaching the 
insulator U Uc2 the quasiparticle weight Z vanishes 
linearly.^'^^'^^ The weight in the first region vansihes like 
Z'^, the weight in the second region like Z while the 
weight in the third region, though small, persists. It is 
also found in the insulator. 

A shoulder occurs in the metallic charge response at 
the same energies as the sharp feature found previously 
at the inner band edges. ^^'^^ The weight in the shoulder 
scales like Z^ so that we are led to the conclusion that 



the shoulder is a consequence rather than the cause of 
the sharp feature in the single-particle propagator. 

The spin susceptibility x^'^"^ in the insulating phase is 
found to be dominated by a strong 6 peak which would 
correspond in finite dimensions to dispersive magnetic 
excitations. In infinite dimensions in a paramagnetic in- 
sulator it happens to be at zero frequency. The peak 
must be seen as a particle-hole bound state. 

Besides this peak only a very weak continuum is found 
at higher frequencies. Hence the localized magnetic mo- 
ment is only very weakly reduced by charge excitations. 
The Mott insulator is governed by very well localized 
spins as long as it exists. 

In the metallic phase, the spin response at higher fre- 
quencies displays again only a continuum of small spec- 
tral weight. Upon increasing interaction the spin spectral 
density is dominated by a pronounced peak at low, but 
finite, frequencies. This peak comprises most of the spec- 
tral weight. It constitutes the precursor of the sharp mag- 
netic mode in the insulator as is evidenced by the coin- 
ciding spectral weights of both features at t/ = Uc2 ■ The 
pronounced metallic peak is the signature of an almost 
bound particle-hole resonance which can be seen as the 
emergent magnetic mode (paramagnon or paratriplon) . 
We expect this mode to persist in finite dimensions as a 
dispersive resonance at low, but finite frequencies. 

This concludes the investigation of the zero tempera- 
ture Mott transition at half-filling in infinite dimensions. 
Further investigations away from half-filling are called 
for. Similarly, it would be very important to verify the 
hypotheses derived here for finite dimensions by future 
calculations. 
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